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1.  INTRODUCTION 


The  a)mbustion  of  nitrogen-containing  species  has  been  of  fundamental  interest  for  over  three  decades 
(Gaydon  1974;  Gaydon  and  Wolfhard  1949;  Wolfhaxd  and  Parker  1955).  Much  of  the  interest  is  directed 
towards  understanding  and  predicting  poUutant  fonnation.  Our  interest  in  nitrogen  combustion  stems  from 
a  desire  to  model  propellant  combustion,  primarily  solid  nitramines  or  nitrocellulose-based  types.  A 
fundamental  undeistanding  of  the  detailed  chemistry  is  a  prerequisite  for  the  successful  modeling  of 
propellant  combustion  which  will  lead  to  optimal  propellant  formulation  and  propulsion  performance.  Of 
foremost  importance  are  molecules  such  as  N2O  and  NOj.  since  they  arc  intermediois  oxidizers  formed 
during  propellant  combustion  (Fifcr  1984;  Schroeder  1985).  Thus,  one  of  the  simplest  flante  chemical 
systans  of  interest  to  both  propellant  chemisiiy  and  pollutam  fonnation  is  the  Hj/N^O  system. 

Several  flame  studies  have  been  performai  on  the  Hj/NjO  system.  Fine  and  Evans  <1963)  measured 
NO  concentrations  and  temperature  in  the  burned  gases  of  flames  using  a  number  of  fuels,  including  H2 
with  N2O.  The  results  were  used  to  infer  qualitative  infonnedon  ubout  reactiwts  responsible  for  NO 
formation  unde*'  *5  wide  variety  of  conditions.  In  a  number  of  otlKr  early  studies  (Dixo«»Lewis>  Sutton, 
and  Williams  1964,  1965a.  1965b;  Duval  arsd  Vaut  Tiggclcn  1967),  flame  spells  and  a  few  species  and 
temperature  profiles  were  also  measured.  More  recently.  Cattolica,  Smooke,  ai>d  Dean  (1982)  perfomted 
a  combined  modeling  and  cxperimeraal  study  for  a  flat,  burner  staWlized,  sicichiomctric,  atmospleric 
[jressurc  Hj/NjO  flstme.  La.ser-iriduced  fluorescence  (LIF)  and  absorption  spectroscopy  were  used  to 
obtain  concentration  profiles  of  OH.  NH.  and  NO.  as  well  as  a  temperature  profile.  However,  the 
measurentems  antld  only  be  performed  in  tlte  uf^ser  portions  of  the  flame  zone.  Subsequent  studies  00 
a  similar  burner  by  Vanderhoff  ct  ai.  (1986)  yielded  ptofiics  for  NO.  O2,  N^.  a»*d  tempemurc  using 
Raman  spectroscopy  and  relative  OH  proi.ics  using  UF  for  several  Hj/NjO  flames  renging  front  lean  to 
stoichiometric  mixture  ratios  ai  amtospheric  pressure. 


Generally,  species  profiles  cgjtaittcd  In  the  atmospheric  pressure  studies  were  poorly  resolved. 
HowTCvcr,  an  iremsasc  in  spatial  resolution  can  be  obtained  u  reduced  pressures  since  the  reaction  rone 
expands.  The  first  low-pressure  work  which  combined  an  experimental  and  theoretical  study  was 
pCffoTttcd  by  BaiaMuiin,  Varxlcoren.  and  Van  Hggelcn  (1977)  on  the  sirurrurc  of  a  lean  flaare 

(3i.4%  Hj)  at  40  torr.  Molcoilar  beam  samj^iitg  with  mass  spcciromctric  detection  (MB/MS)  was 
emf^oyed  to  profile  ail  of  the  stable  and  most  of  the  nsactivo  species.  Using  an  assumed  rttcchantsro, 
consisting  of  approximately  a  dozen  reaaions,  several  blmolecular  rate  constants  were  deduced  and  the 
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importance  of  the  NoO+M  decomposition  -caction  was  discussed.  A  concentration  profile  of  ti'^  N>. 
radical  was  not  reported  and.  as  a  result,  the  NH  chemistr>^  was  excluded  from  the  mechanism.  Recently, 
several  hydrogen  flames  (Kohse-Hdinghaus  et  al.  1988;  Ret^berger  et  al.  1988;  Jeffenes  et  al.  1988) 
supported  by  NjO  or  O2  at  low  pressure  have  been  studied  in  which  the  NH  radical  profiles  have  been 
obtained.  For  one  of  the  H^/NjO  flames,  the  temperature  pnafile  as  well  as  the  OH  atid  NH  radical 
profiles  were  obtained  by  LIF  and  modeled  (Jefferies  et  al.  1988).  A  very  brief  dismtsssod  of  tite 
chemistry  is  given,  but  the  major  emphasis  of  that  work  was  to  determine  the  influence  of  tcmperalu/ 
and  majority  species  mixture  on  quenening  of  the  UF  signals  for  these  radicals. 

In  a  modeling  study,  Coffee  (1986)  developed  a  model  for  promixed,  emc-4imemkstiil  flames  arsi 
calculated  flame  speeds  and  species  profiles  for  a  variety  of  K2/N2O  flames.  The  results  were  compared 
to  a  range  of  previously  published  experimental  data  and  revealed  that  his  model  reproduces  Uie  data 
overall,  but  inaccuracies  still  exist  In  particular,  most  of  the  predicted  flame  speeds  were  slightly  slower 
titan  experiment.  Sensitivity  and  rate  analyses  were  used  to  indicate  that  additional  experimental  data, 
paxiiculariy  on  OH  and  NO.  were  niLcded  to  improve  his  model. 


In  addition  to  the  flame  studies.  itK^rc  have  been  several  shock  tube  (Heiuici  and  Bataer  1969;  IVan 
1976:  Dean.  Steiner,  and  Waitg  1978;  Pamidimukkala  aitd  Skiimer  1982;  Hidaka.  Takunta.  and  Suga 
I9S5a.  1985b)  and  bulb  (Baldwin.  Geiltin,  and  Walker  1973;  Baldwin,  Getiiin,  aiKl  Pisistowx;  1975)  studied 
on  the  Ht/NjO  chemical  system.  In  tlte  .shock  tube  studies,  conditions  were  tailored  primanly  to  highlight 
one  or  a  few  of  Uk  dementafy  reaction  steps  and  to  measure  the  corresponding  rate  coeffidems.  The  bttib 
experiments  examined  the  chemical  mechanism  for  this  system  at  low  temperatures. 


Presented  in  this  report  are  tanperature  and  species  profile  tneasurenenu  for  a  ZO-torr  scoi<d’iOfaetriG 
Hj/N^OAAr  burner  .sEabtlised  flame.  MB/MS  artu  LIF  arc  used  to  measure  a  varidy  of  stable  and  reaciiv’s 
radical  spcdcs.  ITie  tow  pressure  atiaiiKd  en^^s  hi^h  spatial  rcsolutioo  profiles  to  be  obtained  through 
the  cmiro  flame  zone  A  vrider  variety  of  species  profiles  ius  been  obuined  than  fats  pfcvtoudy  appeared 
for  any  study  of  this  chemical  system  due  largely  to  the  use  of  the  MB/MS  sysicni.  Tbs  cxperimearal  data 
is  compared  to  caJeuJated  profiles  genenued  using  the  Sa,ndia  Naional  Uborsiafics  fiarnc  code,  PREMiK 
(Kcc  Cl  al.  1991).  The  mechanism  used  to  generate  Use  model  profiles  is  titc  rcstdl  of  a  detailed 
examination  of  the  literatare  cottceming  gas-phase  nitrogen  combustion.  Rate  and  sertstriviiy  analyses 
performed  reveal  ilte  huricactes  of  the  mechanisnt.  The  main  teatures  of  the  mechanism  are  presented  srvd 


reactions  of  particular  interest  highlighted. 
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2.  EXPERIMENTAL 


A  schematic  of  the  experimental  apparatus  is  shown  in  Figure  1 .  Tlie  H2/N20/Ar  low-pressure  flame 
was  supported  on  a  McKeniia  flat  burner  having  a  6-cm  stainless-steel  fritted  plug.  The  plug  is  encircled 
by  another  sintered  metai  frit  through  which  argon  was  flowed,  thus  forming  a  protective  shroud  that 
minimized  mixing  of  any  recirculating  burnt  gases  in  the  low-pressure  chamber.  The  center  plug  was 
water  cooled  to  maintain  a  constant  temperature  as  measured  by  imbedded  Alumel-Chromel 
thermocouples.  The  burner  was  mounted  at  the  center  of  a  cylindrical  stainless-steel  vacuum  chamber 
(D=25  cm,  H=40  cm)  on  a  high  vacuum  feedthrough  flange  which  was  coupled  to  an  off-axis  rotary 
feedthrough  flange.  Tliis  configuration  allowed  the  burner  to  be  scanned  vertically  with  an  accuracy  and 
reproducibility  of  less  than  50  pm  and  horizontally  with  a  rotational  precision  of  better  than  1°. 

The  reactant  gases  (H2  and  N2O)  and  inert  gas  (Ar)  were  of  commercial  high-purity  grade  and  were 
metered  by  MKS  mass  flow  controllers  (cross-checked  with  a  GCA  Precision  Scientific  wet  test  meter). 
The  gas  flow  rates  employed  were  1.6  and  1.4  L/min  (STP)  for  the  reactant  gases  and  diluent,  respectively, 
corresponding  to  a  mass  flow  rate  of  3.41  x  10“^  g/cmVs.  Pressures  ranging  from  5  to  25  torn  were 
maintained  with  a  D-40  Leybold-Heraeus  vacuum  pump  and  measured  with  a  MKS  270  baratron  (390HA 
pressure  head)  interfaced  to  a  MKS  252  pressure  controller  and  exhaust  throttle  valve.  For  this  study,  the 
flame  was  stabilized  on  the  burner  at  a  pressure  of  20  torr.  Under  these  conditions,  the  flame  was 
characterized  by  a  dark  preheat  zone  and  a  pale  yellow  luminous  zone.  The  dark  or  nonluminous  zone 
is  approximately  5  mm  high  and  occurs  at  the  surface  of  the  burner  where  the  gas  temperature  is  relatively 
low.  The  luminous  zone  is  approximately  5-8  mm  high.  The  thickness  of  the  zones  cou'.d  be  altered  by 
changing  the  relative  proportions  of  the  reactants  and  argon  diluent,  as  well  as  the  pressure. 

Temperature  profiles  of  the  H2/N2O  flawie  were  obtained  with  thermocoup!e.s.  Platinum  and  platinum 
with  10%  ihodium  wires  of  125  diameters  were  spot-welded  to  produce  a  thermocouple  junction  and 
mounted  at  180°  on  a  "V"  shaped  holder  equipped  with  a  spring  to  eliminate  asiy  sag  when  placed  in  the 
flame.  The  thermocouples  were  coated  witli  a  beryllium  oxide  (15%)/yttr{um  oxide  mixture  following  a 
procedu.'e  reported  by  J.  H.  Kent  (1970)  to  avoid  surface  catalytic  effects.  Tlie  temperature  was  measured 
as  tlie  b".mer  traversed  the  distance  from  the  thermocouple  in  the  forward  and  reverse  directions. 
Scanning  tiie  burner  in  both  directions  yielded  temperatures  which  varied  by  less  than  ±5  K.  The  average 
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tesaperatuie  was  tJ«n  conected  for  thermocouple  radi^on  losses  using  ihs  following  equation  (Fristfcm 
and  Westenberg  1965), 


AT  =  T^  .T^=«eada^-T„Vk  .  tD 

where  e  is  dte  cmissivlty  of  tlw  coated  Hwnaocouj^e,  <y  is  tte  SteM  BoltzaiaiTi  constant,  d  is  8!>s  dissie&r 
of  the  coated  junction,  k  is  the  thermal  conductivity  of  the  flame  gases,  and  is  300  K.  The  above 
equation  is  ibr  a  sj^iencai  themiocoupie  junction  with  a  small  enough  diameter  that  its  Reynolds  number 
is  much  less  than  1.  A  measured  diameter  of  190  pm  for  the  coated  therrooctmple  junction  and  an 
emissivit]'  value  of  0.6  reported  by  Pcteisor?  (1981)  were  emjrioyed  in  the  corrccdcn.  The  reported 
emissivity  v.due  lias  an  uncertainty  of  ±12.5%  snd  is  in  good  agreement  with  die  value  of  0.64  repotted 
by  Seery  and  Zabieiski  (1977),  Considering  the  uncertainty  in  the  emissivity  value,  the  mardmum  error 
in  the  corrected  temperature  due  to  this  factor  is  approximately  1%.  The  overall  uncertainty  of  the 
temperature  measurements  is  ±  50  K  in  the  region  of  peak  temperature  where  the  correction  terni  is  largest 
(approximately  200  K)  and  within  20  K  in  the  prslreat  region.  The  sodium  line  revcreal  (SLR)  technique, 
described  in  a  previous  publication  (Bernstein  et  al.  1993),  was  used  to  crosscheck  the  thermocouple 
tempeititute  measurements.  Flame  temperatures  were  measured  in  a  20-ton’  CgH^/OyAr  flame  with  the 
.same  thennocouf^e  and  using  the  SRL  technique.  After  correcting  for  the  thermocouple’s  radiation  losses, 
the  temperature  measurements  were  within  t!te  reported  experimental  uncertainty. 

MB/MS  was  employed  to  monitor  the  stable  and  radical  species.  The  mass  s|wct«5raeter  system 
consists  of  an  Extrel  C50  TQMS  inline  triple  quadruple  mass  filter  with  a  concentric-axis  ionizer  as  shown 
schematically  in  Figure  1.  For  this  study,  the  triple  quadruple  was  operated  as  a  single  quadniple  with 
mass  filtering  by  either  the  first  or  third  quadruple.  Flame  gases  were  sampled  through  a  conical  quartz 
skimmer  with  a  200-pm-diameter  orifice.  Ihe  gases  expanded  supeisonically  into  the  first  differential 
vacuum  chamber  with  a  background  pressure  of  5  x  10"^  ton’.  Tlie  expanding  gases  were  then  collimated 
through  a  second  skimmer  (Beam  Dynamics  Model  2,  2-mm  orifice  diameter)  and  passed  into  the 
ionization  region  of  the  first  quadruple,  maintained  at  a  pressure  of  2  x  10"^  torr.  The  beam  was  then 
modulated  at  200  Hz  with  a  tuning  fork  chopper  and  ionized  prior  to  entering  tlie  first  quadruple 
(ionization  energy  was  17.0  ±0.1  eV  arid  electron  current  maintained  at  0.1  ±  0.01  mA  except  where 
noted).  A  beamstop  was  included  to  determine  if  the  modulated  beam  'was  a  molecular  beam  (modulation 
would  cease  with  beamstop  activated)  or  an  effusive  beam  (no  change  in  modulation  noted  with  beam.stcp 
activated).  The  ion  current  was  measured  with  a  continuous-dynode  electron  multiplier  and  processed  with 


g  lock-in  amplifier  (Stanford  Research  System  Model  SR530).  Phase-sensitive  detection  allowed 
discrminatiori  of  background  gases  and  signal  averaging  to  increase  sensitivity. 

The  ccmcemration  of  stable  species  was  derormined  by  diroa  calibration  at  ambient  temperature  by 
a  {Hocedure  reported  by  Peelers  and  Mahi^  (1973).  The  relationship  between  the  signal  intensity  of  a 
given  species  and  the  partial  pressure  can  be  expressed  as 

Ii  =  SjPi.  (2) 

where  is  the  overall  sensitivity  factor  of  the  apparatus  which  takes  into  account  the  ionization  cross 
section.  It  is  assumed  that  varies  for  all  species  as  a  function  of  {»essure,  temperature,  and  mean 
molecular  weight  of  the  sampled  gas  in  the  same  manner.  Therefore,  the  ratio  of  sensidvidcs  of  any  two 
species,  S^/Sj,  remains  constant  at  every  point  in  the  flame  and  is  equal  to  that  obtained  at  ambient 
temperature.  This  as-sumpdo"  was  verified  by  introducing  a  known  concentradon  of  He  in  the  HjflMjO/Ar 
le  w-pressure  f  ,ne  and  measuring  the  (SHe/S^r)  rado  along  the  flame  zone  and  at  ambient  temperature. 
The  ur.cercainty  in  the  stable  species  concentradon  is  estimated  as  ±  10%. 

The  degree  of  perturbation  inttwuced  by  the  sampler  was  checked  by  monitoring  the  OH  radical  using 
laser-induced  fluoicscence.  The  radical  was  monitored  by  employing  die  Ri(6.5)  and  Ri(7.S)  rotational 
lines  of  its  (i,0)  elect  onic  tra^isition  near  281  nm.  These  rotatiwial  transitions  are  least 

sensitive  to  temperature,  parficularly  in  regions  of  large  temperature  gradient,  as  evidenced  by  Equation  3 
(Eckbreth  1988), 

j2  +  J-(k/hcB,)T,^  =  0  ,  (3) 

where  J  i,s  the  rotational  level;  k,  h,  c,  and  By  are  tlr  usual  spectroscopic  consta^'ts;  and  T^y  is  the  average 
flame  temperature.  The  (1,0)  band  was  chosen  '  'er  the  stronger  (0,0)  band  rrear  306  nm  in  order  to 
minimize  self  absorption.  Overall,  the  profiles  obtamed  by  the  two  techniques  are  in  good  t  ^cement. 
Both  tecimiques  yielded  reliable  species  profile  data  in  the  flame  front  and  burned  zant.  However, 
in  the  preheat  zone,  near  die  burner  .«;urface,  the  quartz  probe  perturbs  tire  flame  slightly.  Similar  findings 
were  reported  in  a  previous  study  (Howard  et  al.  1992). 

The  probe  beam  was  provided  by  an  XeG  excimer-pumpec  dye  laser  (Lumonics  HyperEX 
4(X)/HyperDYE  300)  which  was  frequency  doubled  (hypeffRAK  1000)  to  obtain  tiie  required 
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UV  wavelengths.  The  Hnewidth  for  UV  radiation  is  approximately  0.16  cm“^  (fwhm).  The  probe  beam 
was  collimated  with  irises  to  a  diameter  of  approximately  2.0  mm  and  focused  over  the  center  of  the 
burner  vrith  a  500-mm  focal  length  lens.  The  induced  fluorescence,  corresponding  to  the  (1,1)  transition, 
was  collected  90“  to  the  excitation  laser  beam  and  focused  with  a  3(K)-mm  focal  length  lens  onto  the 
entrance  slit  of  a  Hamamatsu  R-955  {^otomultiplier  equipped  v/ith  a  IS-nm  fwhm  interference  Alter 
catered  at  311  nm  (Corion).  The  resulting  signal  was  then  directed  into  a  boxcar  integrator  (Stanford 
Research  System  2S2)  set  for  a  9-ns  gate  width  in  order  to  minimize  the  effects  of  collisional  quenching 
and/or  energy  trarafer  of  the  rotationaily  excited  states.  The  fluorescence  signal  was  monitored  rai  a 
125-MHz  digital  oscilloscope  (Lecroy  9400).  A  PC-AT  computer  was  employed  for  data  acquisition  and 
aitalysis. 

3.  CODE  AND  COMPUTATIONAL  DETAILS 

The  computations  were  performed  with  the  code  PREMIX  (Ver.  1.4),  developed  at  Sandia  National 
Laboratories  (Kee  et  al.  1991).  This  code  uses  CHEMKIN-II  (Ver.  1.8),  a  library  of  user-friendly, 
chemistry-related  subroutines  and  a  chemical  mechanism  interpreter  code  designed  to  calculate  problems 
involving  elementary  gas  phase  kinetics  (Kee,  Rupley,  and  Miller  1989).  The  equations  for  a  one¬ 
dimensional,  pnemixed,  laminar,  steady-state  flame,  which  includes  complex  chemistry  as  v,cii  as:  gas-phase 
transport  and  thermodynamic  properties,  are  solved  to  predict  species  concentration  temperature 
profiles  as  a  function  of  distance  above  the  burner  surface.  This  code  can  be  used  to  perform  three 
different  types  of  flame  calculations;  1)  free  flame,  2)  burner  stabilized  with  normal  boundary  conditions, 
and  3)  burner  stabilized  with  infinite  H-atom  recombination  rate  to  H2  at  the  burner  surface.  The  first  two 
cases  are  readily  available  in  the  code  while  the  third  case  is  achieved  by  modifying  the  code’s  boundary 
equations  (Wamatz  1978)  for  hydrogen  boundary  conditions  (HBC)  and  has  been  implemented  in  the 
calculations  presented  in  this  repoit.  For  the  present  conditions,  the  solutions  determined  using  either 
HBC  or  normal  boundary  conditions  are  identical  except  for  the  H  atom  profile  close  to  the  burner  surface 
(<5  mm).  The  result  for  normal  boundary  conditions  is,  of  course,  that  the  H  atom  concentration  is  not 
zero  at  the  burner  surface.  The  calculations  also  include  thermal  diffusion  for  light  species  (H  and  H2). 
A  solution  is  regarded  as  acceptable  when  the  following  conditions  are  sadsfled:  1)  the  flame  front  is  far 
enough  from  the  hot  boundary  that  the  boundary  conditions  do  not  affect  the  solution,  2)  the  error 
tolerances  arc  sufficiently  small,  and  3)  increasing  the  number  of  grid  points  will  not  alter  the  solution. 
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PREMIX  extracts  important  gas-phase  thermodynamic  and  transport  properties  for  the  SfKcies  being 
considered  from  two  databases  provided  by  Sandia  National  Laboratories  (Kee,  Rupley^  and  Miller  1987; 
Kee  et  ai.  1988).  The  reactions  employed  in  the  mechanism  are  written  in  the  directiwi  in  which  their  rate 
expressions  are  more  accurately  kirown.  We  define  this  direction  as  the  forward  direction.  At  the  user’s 
discretion,  CHEMKIN  computes  the  rate  constants  of  the  reverse  reactions  by  using  themiodynamic 
properties  and  forward  rate  coefficients.  As  a  result,  thermodynamic  database  can  be  important  (Martin 
and  Brown  1989).  In  the  present  wodc,  several  important  NH  reactions  have  been  reversed.  The 
thermodynamic  properties  for  NH  are  well-known  with  the  exception  of  the  heat  of  formation,  which  has 
been  a  controversial  subject  The  value  used  in  the  database  (Kee,  Rupley,  and  Miller  1987),  however, 
is  very  close  to  the  recommended  value  found  in  a  recent  critical  review  (Anderson  1989).  Version  1.4 
of  ratEMIX  incorporates  the  calculations  of  species  transport  properties  using  the  transport  property 
formulation  (method  VI)  of  Coffee  and  Heimeii  (1981,  1983)  to  obtain  species  transport  effects. 

Recent  work  by  Smith  (1992)  has  demonstrated  that  expansion  of  the  flame  area  above  the  burner  can 
have  a  very  important  effect  on  the  species  profiles,  especially  the  thickness  of  the  flame  zone.  This  is 
due  to  the  residetiice  time  of  the  reacting  gases  (time  it  takes  the  gases  to  travel  a  given  distance)  which 
can  be  strongly  affected  by  expansion.  For  present  conditions,  visual  inspection  of  the  flame  indicated 
that  the  expansion  ratio  is  no  larger  than  1.2  per  cm.  Test  calculations  incorporating  this  value  as  an 
upper  limit  in  the  flame  code  showed  that  the  computed  profiles  are  altered  very  luUe  by  this  effect. 

A  post  processor  was  written  for  the  PREMiX  code  to  analyze  heat  release  as  well  as  sensitivities  and 
rates  of  reactions  at  discrete  distances  above  the  burner  surface.  The  postprocessor  calculates  the  net 
contribution  of  each  reaction  to  the  formation  or  removal  rates  of  given  species  and  sorts  the  reactions  in 
decreasing  order  of  importance.  Tlie  PREMIX  code  calculates  raw  sensitivity  coefficients  for  each 
reaction  and  species,  which  are  nonnalized  by  Uie  post  processor  according  to  the  equation, 

Suc=  A/X,,.m  •  (4) 

where  is  the  normalized  sensitivity  coefficient,  Aj  is  the  Arrhenius  A  coefficient  of  reaction  i,  and  Xj^^^ 
is  the  maximum  mole  fraction  of  species  k.  These  normalized  sensitivity  coefficients  are  sorted  to 
determine  which  reactions  have  the  greatest  effect  on  a  given  species  concentration. 
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The  PREMIX  flame  code  allows  one  to  choose  betwe^  using  the  measured  temperature  profile  as  a 
fixed  input  to  the  problem,  thus  generating  the  species  profiles  for  comparison  with  experiment,  or  to  solve 
the  energy  equation,  thus  obtaining  a  predicted  temperature  profile  as  part  of  the  solution.  The  energy 
equation  in  PREMIX  only  accounts  for  conductive  heat  losses  from  the  flame  to  the  burner  surface.  Both 
a^)roaches  were  tried.  Although  the  fixed  infxit  approach  was  used  to  model  the  experimental  species 
concentration  profiles,  two  important  points  were  found  by  examining  the  results  obtained  using  the  energy 
equation.  First,  a  temperatur?  of  approximately  2,630  K  is  predicted  for  the  burned  gas  just  past  the  flame 
zone,  approximately  200  K  smaller  than  that  predicted  at  a  very  long  distance  fiom  the  burner.  That  is, 
the  postflame  temperature  is  predicted  to  overshoot  its  va'ue  at  equilibrium.  A  temperature  overshoot  of 
similar  magnitude  has  been  observed  by  Smith  (1993)  for  a  low-pressure  H2/N2O  flame.  For  our  flame 
system,  an  overshoot  is  also  predicted  for  the  corresponding  free  flame  calculation  where  the  correct 
adiabatic  flame  temperature  is  stained  at  large  distances.  These  unusual  temperature  overshoots  are  due 
to  the  predicted  concentrations  of  H,  O,  and  OH  radicals  in  the  flame  zone.  The  typical  radical  overshoot 
is  not  predicted  for  these  flames.  Instead,  the  radical  concentration  is  lower  than  the  equilibrium  value 
and  rises  slowly  in  the  postflame  region.  The  second  point  in  solving  the  energy  equation  is  that  it  results 
in  a  predicted  final  flame  temperature  that  is  approximately  50Q-6G0  K  larger  than  that  measured 
experimentally.  Not  surprisingly,  the  predicted  flame  zone  width  is  much  smaller.  A  400-500  K 
difference  is  also  observed  by  Smitli  (1993)  in  an  independent  study  on  tlie  H2/N2O  flame.  Large 
differences  are  predicted  for  several  other  fuel/oxidizers  as  well.  These  results  suggest  that  nonnegligiblc 
gas  phase  heat  losses  (for  example,  radiative)  to  the  surroundings  are  not  properly  accounted  for  in  the 
energy  equation.  Since  the  predicted  temperature  is  so  different  from  experiment,  this  approach  cannot 
reasonably  be  used  to  model  tlie  results.  It  is  clear  that  further  study  of  tlie  heat  loss  mechanism(s)  is 
needed.  For  our  system,  it  should  be  noted,  however,  that  the  magnitudes  and  relative  ordering  of 
sensitivities  of  the  species  profiles  to  reaction  rate  coefficients  were  found  to  be  similar  for  both  methods. 
This  means  that,  peitiaps  suiprisingiy,  the  computed  species  profiles  arc  nearly  equally  sensitive  to  Uie 
kinetics  parameters  for  either  approach. 

4.  RESULTS  AND  DISCUSSION 

4.1  Chemical  Mechanism.  The  gross  chemical  stnicture  of  H2/N2O  flames  can  be  characterized  by 
a  number  of  important  reactiotu  wtiich  have  been  discussed  in  previous  studies  (Cattolica,  Smookc,  and 
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Dean  1982;  Balakhnin,  Vandoonsn,  and  Van  Tiggelen  1977;  Co^ee  1986).  The  key  r^cal-pnsducing  siep 
promoting  the  combustion  under  most  conditions  is: 


N20  +  M  =  N2  +  0  +  M  .  (5) 

Most  of  the  reactants,  H2  and  N2O.  are  convened  to  products,  H2O  and  N2,  by  die  propaption  stqis: 


NjO  +  H  =  OH  +  N2 


(6) 


H2  +  OHSH2O  +  H  . 


(7) 


A  small  iiaciion  of  the  N2O  is  converted  to  NO,  which  is  observed  at  a  few  mole  percent  In  the  burned 
gases  under  most  conditions.  In  prior  mechanisms,  the  reaction  N2O  +  H  »  NH  +  NO  was  believed  to 
be  responsible  for  most  of  the  NO  production  under  rich  or  stoichiometric  conditions.  Towards  lean 
conditions,  the  reaction  N2O  -t-  O  =  NO  +  NO  assumes  greater  significance  for  the  production  of  NO. 

For  our  study,  a  comprehensive  chemical  mechanism  containing  over  2(X)  reactions  and  more  than  20 
species  was  initially  used  to  model  the  detailed  chemistry  of  the  flame.  Most  of  the  reactions 

and  rate  expressions  were  obtained  from  the  Miller  and  Bowman  mechanism  (Miller  and  Bowman  1989) 
for  gas-phase  combustion  involving  nitrogen  compounds.  Added  to  this  reaction  set  wete  the 
recombination  reactions  of  NH^  from  ammonia  chemistry  (Dean,  Chou,  and  Stem  1984)  and  combustion 
reactions  involving  NO2  and  N2O  from  a  recent  critical  review  of  the  chemical  kinetic  database  for 
propellant  combustion  (Tsang  and  Herron  1991).  A  few  reactions  for  predicting  NO^  formation  and 
destruction  were  obtained  from  a  paper  on  HCN/NO2  combustion  (Thorne  and  Melius  1989).  A  sensitivity 
analysis  perfonned  on  this  comprehensive  mechanism  revealed  reactions  which  were  not  important  for  tiu: 
present  experimental  conditions.  Most  of  these  reactions  were  eliminated,  and  as  a  result,  species  sucti 
£s  H2O2,  NH2,  NH3,  NjH^,  NO2,  NO3,  N2O4,  HONO,  and  HNO3  were  removed  from  the  reaction  set. 
The  resulting  mechanism,  listed  in  Tabic  1.  contains  38  reactions  and  14  species  and  is  med 
the  modeling  results  presented  in  the  following  sections.  WIten  two  reactions  involving  NNH,  whose 
importance  is  discussed  and  discounted  in  the  following  sections,  are  removed,  tiie  resulting  predictions 
of  the  comprehensive  mechanism  are  virtually  identical  to  those  from  our  smaller  meciianlsm. 
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Table  I.  Reaction  Mechanism  Rate  Coefficients  in  the  Fbmi  k  = 


Reaction 

A* 

B 

C 

Reference 

1. 

Hi  +  Oi  =  OH  +  OH 

1.70E13 

0.00 

47,780.0 

Miller  and  Bowman  1989 

2. 

Hi  +  OH  =  HjO  +  H 

2.16E08 

U1 

3,430.0 

Michael  ami  Suthcriand  1988 

3. 

H,  +  0  »  H  +  OH 

5.06E04 

2.67 

6290.0 

Miller  and  Bowman  1989 

4. 

H  +  H  +  M  =  H,+  M 

HyOy  HjO/O/ 

1.00E18 

-1.00 

0.0 

Miller  and  Bowman  1989 

4a, 

H  +  H  +  Hi  =  Hi  +  Hi 

9.20E16 

-0.60 

0.0 

Miller  and  Bowman  1989 

4b. 

H  +  H  +  HjOsHi+HiO 

6.00E19 

-125 

0.0 

Miller  and  Bowman  1989 

5. 

H  +  Oi  =  OH  +  0 

3,52E16 

-0.70 

17,070.0 

Masten  et  al.  1990 

6. 

0i  +  H  +  M  =  H02  +  M 
(V1.0/  N,fl.3/  HiO/18.6/  Hj/2.9/’ 

3.61E17 

-0.72 

0.0 

Miller  and  Bowman  1989 

7. 

0  +  0  +  M  =  0,  +  M 

1.89E13 

0.00 

-1,788.0 

Miller  and  Bowman  1989 

8. 

HOj  +  OH  =  H,0  +  Oi 

7.50Ei2 

0.00 

0.0 

Miller  and  Bowmai  1989 

9. 

HOi  +  H  *  OH  +  OH 

1.69E14 

0.00 

874.0 

Tsang  a?ai  Hamj^oo  5986 

10. 

HOi  +  H  =  Hi  +  Oj 

6.63E13 

0.00 

2,126.0 

Tsang  and  Ha'.n?;aon  1986 

11. 

HO,  T  0  »  OH  +  0, 

1.40E13 

0.00 

1,073.0 

Miller  /ind  Bowman  1989 

12. 

OH  +  OH  »  HjO  +  0 

6.00E08 

1.30 

0.0 

Miller  and  Bowman  1989 

13. 

H  +  OH  +  M  a  H,0  +  M 

H,0/5/ 

1.60E22 

-2.00 

0.0 

Miller  and  Bowman  1989 

14. 

H  +  0  +  M  =  0H  +  M 

H,0/5/ 

6,20EI6 

-0.60 

0.0 

Miller  and  Bowman  1989 

15. 

N,0  +  H  a  OH  +  N, 

2.53E10 

0.0 

4,550.0 

Marshall  ct  al.  1989 

N,0  +  H  =.  OH  +  N, 

2J13E14 

0.0 

16,750.0 

Marshall  ct  al.  1989* 

16. 

N,0  +  M  »  N,  >  0  +  M  7.20E17  >-0.73 

No/oy  N,o/oy  hjoj  ojoj  hjo/oj  njoj  m  a)/ 

62,800.0 

Tsaixg  and  HetTon  1991'* 

i6a. 

N,0  +  NO  =>  N,  0  +  NO 

7.20E17 

-0.73 

62,800.0 

Tsang  and  Herron  1991* 

I6b. 

NjO  H,  ^  N,  +  0  +  H, 

7.20E17 

-0.73 

62.800.0 

Tsang  and  Henon  1991* 

16c. 

N,0  +  0,  oNj  +  O  +  O, 

7,20E17 

-0.73 

62300.0 

Tsang  and  Heiroo  1991* 

*  Unitti  A»CTn-mal-s«c-K,  Cacil/mole.  Lovv-presture  limits  used  for  unimolocuUr  docomposiUon  and  rocombiiulioci. 

*  Third  body  effidcocics. 

°  Tha  reaction  rit«  coefficient  is  computed  by  the  sum  of  the  two  cxpresrtcns. 

TTiird  body  efficiencies  for  Ar,  N^O,  N,  Grom  Bsukh.  Drystkle.  Ktd  Ifome  (1973)  and  H]0,  H,  ind  NO  froni  best 
estimate. 
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Table  1.  Reaction  Mechanism  Rate  Coefficients  in  the  Fonn  k  =  (continued) 


Reaction 

A 

B 

c 

Reference 

16d. 

NjO  +  Nj  =  N2  0  +  Nj 

7.20E17 

-0.73 

62,800.0 

Tsang  and  Herron  1991“ 

16c. 

NjO  +  NjO  =  Nj  +  0  +  NjO 

3.60E18 

-0.73 

62,800.0 

Tsang  and  Herron  1991* 

I6f. 

N2O  +  HjO  =  N2  +  0  +  H2O 

3.60E18 

-0.73 

62300.0 

Tsang  and  Herron  1991“ 

16g. 

N2O  +  At  =  N2  +  0  +Ar 

4.82E17 

-0.73 

62,800.0 

Tsang  and  Herron  1991“ 

17. 

N2O  +  0  =  N2  +  O2 

1.00E14 

0.00 

28,000.0 

Tsang  atid  Herron  1991 

18. 

N2O  +  0  =  NO  +  NO 

6.60E13 

0.00 

26,600.0 

Tsang  and  Herron  1991 

19. 

NjO  +  OH  =  HO2  +  N2 

2.00E12 

0.00 

10,000.0 

Miller  and  Bowman  1989 

20. 

NH  +  NO  =  N2O  +  H 

2.94E14 

-0.40 

0.0 

Miller  and  Melius,  in  press** 

NH  +  NO  =  N2O  +  H 

-2.16E1 

-0.23 

0.0 

Miller  and  Melius,  in  press'* 

21. 

NH  +  NO  =  Nj  +  OH 

2.i6E13 

-0.23 

0.0 

Miller  and  Melius,  in  press 

22. 

NH  +  O2  =»  NO  +  OH 

7.60E10 

0.00 

1330.0 

Mertens  et  al.  1991 

23. 

NH  +  Oj  =>  HNO  4-  0 

3.89E13 

0.00 

17,885.0 

Mertens  ct  al.  1991 

24. 

NH  +  OH  =*  HNO  +  H 

2.00E13 

0.00 

0.0 

Miller  and  Bowman  1989 

25. 

NH  +  OH  =  N  +  HjO 

S.OOEll 

0.30 

2,000.0 

Miller  and  Bowman  1989 

26. 

NH  +  N  =  N2  +  H 

3.00E13 

0.00 

0.0 

Miller  and  Bowman  1989 

27. 

N  +  Hj  »  NH  +  H 

1.60B14 

0.00 

25,140.0 

Davidson  attd  Hanson  1990 

28. 

NH  4.  0  »  NO  +  H 

5.50E13 

0.00 

0.0 

Mertens  cl  al.  1991 

29. 

NH  4-  0  =  N  +  OH 

3.72E13 

0.00 

0.0 

Mertens  ct  al.  1991 

30. 

NH  +  NH  a  Nj  4  H  +  H 

5.IOE13 

0.00 

0.0 

Mertens  et  al.  1989 

31. 

NH  +  M  =  N4-H4-M 

2.65E14 

0.00 

75314.0 

Mertens  cl  al.  1989 

32. 

HNO  4-  OH  =  NO  +  HjO 

4.80E13 

0.00 

990.0 

Tsang  and  Henron  1991 

33. 

NO  4-  N  «  Nj  4-  0 

3.27S12 

OJO 

0.0 

Miller  and  Bowman  1989 

34. 

N0  4-M  =  N4-04-M 

1.40E15 

0.00 

148.430.0 

Tsang  and  Herron  1991 

35. 

NO  4-  H  4-  M  »  HNO  4-  M 

9.00E19 

-1.30 

735.0 

Tsang  and  Herron  1991 

36. 

NO  4-  H  «  N  4-  OH 

1.68E14 

0.00 

47370.0 

Tsang  and  Henon  1991 

37. 

NO  +  0  =  N  +  Oj 

3.80E09 

1.00 

41375.0 

Tsang  and  Herron  1991*^ 

38. 

N  4  NjO  a  Nj  4-  NO 

1.00E13 

0.00 

19370.0 

Hanson  snd  Salimion  1985 

*  Third  body  efficiencies  for  Oj.  Ar,  NjO.  Nj  froni  Btuich,  OryttUJc,  end  Honw  (1973)  tnd  H2O.  *l>d  KO  froHi  best 
estimtic. 

**  The  rc«aion  rite  coefficiem  is  computed  by  Lhc  sum  cf  the  iv?o  Mprewic-^. 

^  The  exprcsiion  shovm  is  from  Hinson  *nd  Silimian  (J985X  »?hieh  is  recommended  in  Tting  ir«d  Herron  (1991).  There  is 
t  trinscripoon  cnor  in  Tsing  and  Herron  (1991). 
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The  mechanism  was  reviewed,  and  reactions  which  arc  important  for  the  present  condititms  were 
carefully  examined.  The  rate  expression  reported  by  Marshall  and  coworkers  (Marshall,  Fontijn,  and 
Melius  1987;  Marsha,  Ko,  and  Fontijn  1989)  for  the  reaction  NjO  +  H  =  OH  +  Nj  (R15)  was 
incorporated  in  the  mechanism.  Those  experimental  and  theoretical  studies  demonstrated  that  the  reaction 
rate  expression  possesses  an  upward  curvature  in  the  Arrhenius  plot  at  low  temf^rature  due  to  quantum 
mechanical  tunneling.  Therefore,  extrapolating  previously  reported  high  temperature  expressions  to  low 
temperatures  could  produce  inaccurate  conclusions.  Rate  expressions  resulting  from  the  recent  work  of 
Hanson  and  coworkers  (Meitens  et  al.  1991;  Da’ddson  and  Hanson  19W;  Mertens  et  al.  1989)  for  a 
number  of  N,  NH,  and  NO  reactions  (R22,  R23,  R27-R31)  studied  at  high  temperatures  have  also  been 
incorporated.  In  particular,  R27  and  R33  have  also  b^n  examined  by  Koshi  et  al.  (1990),  using  similar 
techiuques.  However,  Davidson  and  Hanson  (1990)  reported  careful  calibration  of  experimental 
diagnostics,  and  their  error  limits  are  much  smaller  than  those  reported  by  Koshi  et  al.  (1990).  Therefore, 
the  rate  expression  repotted  by  Hanson  and  coworkere  was  chosen  for  R27  and  is  extrapolated  to  lower 
temperatures  for  the  pns^sent  work.  The  expression  of  Miller  and  Bowman  was  retained  for  R33  beermse 
it  is  in  good  agreement  with  both  the  higli  temperature  data  (Davidson  and  Hanson  1990)  and  results  of 
low  temperature  critical  reviews  (DeMore  et  al.  1987;  Atkinson  et  al.  1989).  Other  reactions  important 
for  our  conditions  will  be  discussed  in  Uic  following  sections. 

^■2  CQtnpari.son  Between  Model  atKi  Extierimem.  Presented  in  Figure  2  are  Uie  experimentally 
measured  temperatures  of  a  stoichiometric,  20-iorT  Hj/N^O/Ar  fianre.  The  data  arc  fitted  by  a  sigmoidal 
function  (Miher  and  Kotlar  1986)  and  tire  resiling  profile,  together  with  an  exirapolatcid  burner  surfaa: 
temperature  of  610  K.  is  used  as  a  fixed  input  for  the  fianic  code.  ITtc  computed  profiles  for  N,0.  H,. 
Nj,  K2O.  aitd  Ar  gerterated  using  dte  mechanism  with  hydiogcn  boui'idary  coruliuons  are  cooipaied  to  the 
MB/MS  experimental  pmfiles  iii  Figure  3.  Tltc  flame  frocu  e.xxciuJs  from  Ure  bunocr  surface  to 
approximately  15  mm  above  the  Iximer,  as  seen  by  the  disappearance  of  die  teacianis  N^O  and  Hj.  In 
the  burned  gas  region,  dre  .'.ompuced  arrd  measured  profiles  for  N2O  deviate  slighUy,  white  the  Kj  profiles 
are  in  agreement.  However,  near  the  burner  surface,  the  NjO  profiles  compare  wv^ll.  but  the  computed 
conceiuraiions  for  H2  are  much  lower  Uian  die  experimerual  values.  The  quart?,  skimmer  may  iniiibii 
diftusioir  near  the  burner  surface  which  would  explain  die  higher  cxpciimcntal  H2  (Xjuccsuration  value. 
The  computed  profiles  for  the  products  Nj  and  H^O  agree  with  the  experimental  results  over  tk;  entire 
flame  zone,  llic  model  predicts  the  absolute  concmraiion  values  well  at  both  Uie  burner  surface  and  in 
the  burned  gas  regiw.  However,  the  model  prediction  for  HjO  in  the  burned  gas  region  is  slightly  larger 
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TEMPERATURE  (K) 


than  the  experimental  value  which  may  result  from  difficulties  in  the  calibration  of  Lite  MB/MS  H^O 
signal.  Overall,  the  agrccmcnt  between  the  measured  and  predicted  majority  species  profiles  is  acceptable. 

Presented  in  Fipre  4  are  the  calculated  and  experimental  concentrallcn  profiles  for  NO  ami  NH.  The 
absolum  value  of  the  predicted  NO  concentration  agrees  well  witli  the  experimental  resula  in  the  peetflsme 
tvsgion,  even  though  there  may  be  a  slew  reaction  consuming  NO  in  this  region  mining  in  the  model  The 
NH  profile  i$  normalized  to  the  model  cotKentration  profile  since  absolute  mea»ur^e»is  wem  not 
obtained  for  NH.  Althou^  the  model  and  cxt^'imental  shapes  of  the  NH  profiles  are  similar,  the 
modeled  NH  profile  is  shiiled  approxima^ly  2  mm  further  away  from  the  burner  surface  than  the 
cxpefimentally  measured  one.  It  is  difficult  to  cetermine  whetlter  tlis  is  tits  result  of  a  slight  perturhfdion 
induced  by  the  MB/MS  technique  or  inaccuracies  in  the  mechanism. 

Experimental  profiles  measured  using  botli  MB/MS  and  UF  for  the  OH  radical  are  shown  in  Figuro  5 
along  with  the  model  results.  Absolute  values  were  not  detennined  in  either  expcriiiacnt.  ihercforo,  the 
experimental  profiles  arc  normalized  to  ?he  calculfed  cotteentmkMis.  Tlw  relative  profiles  for  ixk,'  UF 
experimental  results  agree  very  well  with  ilte  mtxk!.  The  profile  for  the  MB/MS  agr^  reasonably  witlr 
the  model  profile  except  near  the  burner  surface  (1  to  6  mm)  wl®re  the  MB/MS  shows  a  reduced 
concentration  at  4  ram.  Ihis  feature  is  probably  m  anifset  of  dte  MB/MS  technique  since  die  UF  profiles 
do  not  show  Uiis  arKwnaly.  The  UF  tcclmique  is  less  inuusivc  and  care  has  been  taken  to  avx^  possible 
perturbing  effects  such  as  saturation,  temperature.  varifiUotr.  tmd  quenching. 

Figure  6  shows  the  csperimcmally  measured  profiles  of  th:  radica?  spedos  0  and  H  ivocmaliz!^  to 
model  profiles.  The  0-atora  profile  shown  is  an  average  of  two  scans  tccenkd  lutii:^  an  ckorotr  ssis^gy 
of  15.3  ±  0.1  cV.  just  below  tire  tlueshold  for  forming  O*  from  NjO  (Collin  aeid  Lossieg  1958).  Ovcnsil, 
the  agrecntcni  between  esperimem  arrd  model  profiles  is  good  coasklerirtg  the  rather  poor  sifaRl-u>*nDisc 
in  the  experimeni.  In  the  case  of  the  K-aiotn  profiles,  agrocmem  betwesen  ihc  eoqrerimcni  >iad  model  is 
fair. 


Presented  in  Figure  7  are  the  calculated  profiles  of  Oj  srrd  HOj.  Various  rate  coeffidcnts  for  nstciioo 
R19  were  used  .since  the  exact  expoKsion  is  not  w^l  cstaWishod,  The  ertperimemai  profile  for  Oj  is 
predicted  reasonably  wicll  for  two  of  the  model  profiles  a.nd  will  be  discussed  later.  Tite  KOj  r^cal. 


CONCENTRATION  (xl 


CONCENTRATION  (xlO  mole/cc) 


7.  MsasHred  pd  competed  profiles  of  O^;  and  computed  profiles  of  HQ.,.  Ttw  4sn|iit.. 
■SirO^  was  me^tsured.  The  rate  coefficient  for  N,0  +  OH  =  HO^  +  N.  (R19^  was  ohanpM 
fignUiiejnechpism  (pljd  line)  to  dtc  critical  review*s  (Tsane  and  Heiron  1991)  unaer 
HmU  (dotted  line)  pd  zero,  omitted  from  the  mechanism  fdashed  HneV 
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however,  could  not  be  detected  under  our  experimental  conditions.  Unfortunately,  firm  conclusions 
regarding  the  chemistry  cannot  be  drawn  since  the  lower  limit  of  detection  for  this  species  is  not  known. 

The  effects  of  temperature  on  the  predicted  species  concentration  profiles  were  also  examined  since 
the  experimental  temperature  data  was  used  as  a  fixed  input  to  the  fiame  code  calculations.  A  ±5% 
difference  in  the  temperature  profile,  twice  the  experimental  uncertainty,  results  in  a  ±10-15%  difference 
in  the  computed  NO  and  O2  burned  gas  concentrations,  a  ±20-30%  difference  in  the  H,  O,  and  OH 
concentrations,  and  approximately  ±5%  difference  in  the  NH,  N,  HNO,  and  HOj  trace  intermediates.  The 
flame  zone  width,  as  determined  from  the  computed  species  profiles,  changes  only  slightly  compared  to 
the  observed  width,  approximately  ±0.5-1  mm.  Thus,  only  the  computed  H,  0,  and  OH  profiles  have 
much  dependence  on  temperature.  The  significance  of  these  changes  cannot  be  examined  since  the 
absolute  concentrations  of  these  species  were  not  measured.  Even  the  20-30%  changes  are  not  large 
compared  to  typical  error  limits  in  measured  absolute  radical  concentrations. 

4.3  Flame  Stnieture  Analysis. 

4.3.1  Overview  of  the  Nitrogen  Chemistry.  Analysis  of  the  postprtjcessor  results  yields  reaction 
pathway  diagrams  constructed  using  the  mechanism  in  Table  1.  The  reaction  pathway  diagram  depicting 
the  nitrogen  chemistry  occurring  at  approximately  7  mm  above  the  burner  surface  is  shown  in  Figure  8. 
This  point  is  halfway  through  the  fiame  zone  and  is  calculated  to  be  the  position  of  maximum  chemical 
heat  release  in  the  flame.  Qualitatively  this  diagram  is  not  expected  to  change  much  at  various  heights 
in  the  flame  zone  due  to  the  single  stage  nature  of  the  flame.  The  numbers  in  paientheses  are  the  relative 
rates  of  the  various  reactions,  nonnalized  to  KK)  for  the  reaction  N2O  +  H  =  OH  +  Nj  (R15)  representing 
the  fastest  step  consuming  N2O.  The  initiation  step,  N2O  +  M  =  N2  +  O  +  M  (R16),  and  the  reaction 
N2O  +  OH  =  HO2  +  N2  (R19)  are  also  important  reactions  in  the  direct  conversion  of  N2O  ro  N2. 
According  to  the  diagram,  prediction  of  the  concentration  of  NO  is  quite  complex  because  a  portion  of 
the  NH  foraied  by  N2O  +  H  =  NH  +  NO  (-R20)  reacts  witli  H  or  OH  to  form  N  atoms,  which  in  turn 
primarily  react  with  NO  to  form  N2.  However,  there  are  a  number  of  other  competing  pathways  involving 
NH  and  N  atoms  which  may  result  in  the  formation  of  either  additional  NO  or  N2.  Surprisingly,  most 
of  the  rate  coefficients  for  the  reactions  which  affect  the  computed  NO  concentration  have  been  carefully 
studied. 
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Infomiation  concerning  ?lic  overall  combustion  rate  can  be  found  by  examining  sensitivities  of  a  major 
species.  For  example,  the  sensitivity  plot  for  N2O  in  Figure  9  was  calculated  using  the  reaction  set  listed 
in  Table  1  and,  as  expected,  the  radical-producing  reaction  R16  has  the  largest  effect  on  the  N^O 
concentration.  Sensitivities  of  the  propagating  steps  ^12  and  R15  are  also  ’arge. 

As  mentioned  previously,  a  sensitivity  analysis  was  also  run  on  the  comprehensive  i:iechanism 
(>200  reactions).  Surprisingly,  the  overall  combustion  rate  was  found  to  be  highly  sensitive  to  the 
reactions  N2O  +  H  =  NNH  +  O  (R39)  and  NNH  +  M  =  N2  +  H  +  M  (R40).  The  overall  combustion  rate 
was  sensitive  to  these  two  reactions  since  the  sequence  is  chain  branching  and  produces  two  radicals.  The 
sensitivity  coefficient  for  R39  had  approximately  the  same  magnitude  as  that  for  the  initiation  step, 
N2O  +  M  (R16).  Additionally,  the  rate  of  radical  pool  buildup  by  these  two  paths  was  predicted  to  be 
nearly  equal.  The  rate  expression  for  -R39,  taken  from  the  Miller  and  Bowman  mechanism  (Miller  and 
Bowman  1989),  is  an  upper  limit  estimate  (Miller  1993)  (l.(X)  x  10^'^cmVmoI-s).  R40  was  taken  from 
Miller  et  al.  (1983)  and  collider  efficiencies  for  R40  were  assumed  to  be  identical  to  R6.  V/lien  these 
reaction?  were  included  in  the  medianism  the  predicted  flame  speed  increased  and  the  predicted  shapes 
of  the  species  profiles  are  compressed  slightly  (10-15%)  towards  the  burner.  The  absolute  concentrations 
in  the  burned  gas  region  for  all  species  change  only  slightly  (<5-10%)  except  for  lire  H,  0,  and  OH 
concentrations  which  increase  by  ~30-75%.  Unfortunately,  it  is  not  possible  to  determine  whether  these 
reactions  are  important  based  on  the  present  experimental  profiles  since  botli  sets  of  predictions,  with  or 
without  R39  and  R40,  agree  within  error  limits  of  the  experimental  results. 

For  the  reaction  of  H  with  N2O,  Marshall  and  coworkers  (Marshall,  Fontijn,  and  Melius  1987; 
Marshall,  Ko,  and  Fontijn  1989)  present  upper  limit  estimates  based  on  thermodynamic  arguments  for  the 
rate  constant  expression  for  the  pixxlucts  NH  +  NO  (-R20)  and  for  the  branching  ratio  to  NNH  +  O  (R39) 
at  2,000  K.  The  activation  energy  for  R39  is  taken  to  be  no  smaller  than  the  endothcrmicity  of  tlie 
channel.  Following  their  methods,  the  upper  limit  rate  expression  (k  =  5  x  10''*e^~^‘®^^\mVtnol-s) 
for  the  NNH  +  0  channel  is  approximately  a  factor  of  30  lower  than  what  is  obtained  by  reversing  the 
expression  of  Miller  and  Bowman  in  the  temperature  range  of  interest  Wlicn  this  smaller  rate  coefficient 
is  used  in  the  mechanism,  tire  reaction  has  little  significance,  even  for  higlier  temperature  conditions  (c.g., 
(j)  =:  1.0,  atmospheric  pressure  free  flame,  no  diluent).  ITie  apparent  inconsistency  upon  reversing  rate 
coefficients  is  possibly  due,  in  part,  to  uncertain  Uicrmodynamic  data  for  NNH.  Therefore,  rcaciitms  R39 
and  R40  were  excluded  from  our  mechanism. 
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N^O  NORMALIZED  SENSITIVITIES 
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4.3.2  Conversion  of  N^O  to  Nj.  As  discussed  earlier,  the  reaction  NjO  +  H  =  OH  +  Nj  (R15)  plays 
the  prominent  role  in  the  direct  conversion  of  NjO  to  Nj.  The  remainder  of  the  direct  conversion  results 
primarily  from  reactions  NjO  +  OH  =  HO2  +  Nj  (R19)  and  NjO  +  M  =  Nj  +  0  +  M  (R16).  It  has  been 
suggested  by  a  number  of  researchers  (Miller  and  Bowman  1989;  Schoffield,  Vandooren,  and  Van 
Tiggelen  1986;  Tsang  and  Herron  1991;  Yetter  et  ai.  1991)  that  R19  may  be  important  in  several  systems; 
however,  its  rate  constant  is  poorly  established.  Yetter  et  al.  (1991)  found  it  necessary  to  include  this 
reaction  to  model  their  Hj/NjO  flow  system  experiments  at  intermediate  temperatures  (925-1,073  K). 
Under  these  conditions,  the  effects  of  this  reaction  on  the  overall  combustion  are  inhibitive,  presumably 
because  HOj  is  a  slowly  reacting  radical  compared  to  H,  0,  and  OH.  and  also  because  formation  of  HOj 
results  in  removal  of  radicals  from  the  system  by  HOj  +  H  =  Hj  +  Oj  (RIO).  Use  of  the  upper  limit  rate 
expression  for  R19  from  the  critical  review  of  Tsang  and  Herron  (1991)  resulted  in  a  strong  inhibiting 
effect  which  did  not  properly  model  the  flow  system  results.  Reasonable  agreement  was  found  when  using 
the  expression  of  Miller  and  Bowman  (1989).  a  factor  of  3-4  smaller.  Miller  and  Bowman  derived  their 
expression  to  help  explain  concenuation  profiles  measured  in  an  NHyOj  flame  by  Bian,  Vandooten.  and 
Van  Tiggelen  (1986).  Using  their  mechanism.  Miller  and  Bowman  predicted  that  NjO  in  that  flame  is 
formed  primarily  by  NH  +  NO  =  N2O  +  H  (R20)  and  consumed  by  R15  and  R19.  At  that  time,  the  rate 
constant  of  R20  at  high  temperature  was  pwrly  established.  Miller  and  Bowman  stated  lltai  if  subsequent 
work  resulted  in  an  irtcrcase  of  this  rate  constant,  one  possible  way  to  compensate  for  the  increase  in  N2O 
concentration  without  aflecting  seriously  any  of  the  other  predicted  species  profiles  would  be  to  increase 
the  rate  coclTicicnl  of  R19.  As  will  be  discussed  later,  a  slight  increase  at  tlic  high  temperature  end  of 
the  rate  expression  for  R20  is  herein  recommended.  Unfortunately,  this  creates  a  dilemma  because  die 
present  experimental  results  would  be  best  explained  if  the  rate  coefficiciu  for  R19  were  smaller  than  that 
from  the  Miller  and  Bowman  expression. 

The  mtc  expression  listed  in  Table  1  for  R 19  is  from  Miller  and  Bowman  (1989),  Examining  the 
sensitivity  diagram  for  O2.  sfiown  in  Figure  10.  reveals  dial  the  Oj  conenuration  is  sensitive  to  RI9  and 
to  the  branching  ratio  between  R9  and  RIO  since  RIO  forois  O2  directly.  All  of  tire  measured  species 
exhibit  liulc  or  no  .sensitivity  to  R19  except  for  O^.  Figure  7  stows  calculated  Oj  and  HOj  profiles 
employing  different  rate  expressions  for  R19.  The  solid  curve  is  generated  using  Uic  mechanism  in 
Table  1 .  while  tire  other  two  curves  were  generated  by  changing  the  rale  expression  to:  1 )  upper  limit  as 
recommended  by  Tsang  and  Henon  (1991);  and  2)  ?.cro  (tiiat  is.  the  reaction  is  not  included).  If  the  upper 
limit  rate  expression  is  used,  then  the  O2  sensitivity  to  this  rcaciion  itweases,  as  cxp«:-tcd.  Also.  N2O 
sensitivity  to  R9  and  RIO  bccome.s  imporu.ni,  and  inhibiting  effects  arc  rroticcaWe.  As  seen  in  Figure  7, 
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a  large  peak  at  10  mm  is  noticed  in  the  Oj  profile  when  the  upper  limit  rate  expression  is  used.  Tlie 
experimental  profile  obtained  with  the  MB/MS  clearly  does  not  show  a  peak.  Adequate  agreement  between 
the  model  and  experiment  in  this  region  of  the  flame  can  only  occur  while  still  maintaining  the  large  rate 
coefficient  for  R19  if  several  of  the  reactions  to  which  the  Oj  profile  is  most  sensitive  are  changed 
sufuUtaneously  to  the  limits  of  their  uncertainties.  Such  a  drastic  alteration  of  the  mechanism  seems 
unreasonable  even  though  the  other  calculated  profiles  are  only  slightly  affected.  If  reaction  RI9  is 
excluded  from  the  mechanism,  then  the  predicted  and  experimental  O2  profiles  agree  well.  The  above 
results  indicate  that  the  Miller  and  Bowman  rate  expression  for  R](9  is  more  consistent  with  our 
experimental  conditions  titan  the  higher  critical  review  value.  Moreover,  leaving  the  reaction  out  or 
lowering  the  rate  constant  further  yields  results  which  are  more  consistent  with  our  experimental  data. 
However,  without  any  justification  for  lowering  the  e.xpcession  for  R19,  our  recommendation  is  to  retain 
it  in  the  mechanism.  Further  in\  estigaiion  of  this  reaction  is  clearly  needed. 

The  NjO  concentration  and,  hence,  overall  combustion,  is  highly  sensitive  to  the  initiation  aep, 
NiO  +  M  =  Nt  0  +  M  (R16).  as  seen  iri  Figure  9.  To  further  examine  this  reaction,  the  effects  of 
individual  collision  partners  were  included  in  the  reaction  mechanism  (R16a-R16g)  in  order  to  assess 
which  colliders  are  mo^  important,  lire  third  body  efficiencies  for  O2.  .Ar.  N^O.  and  N2  are  fnjra  a 
critical  review  (Baulch,  Drysdale.  ai¥l  Home  1973)  and  the  values  for  HjO.  aiKl  NO  are  csliraaied. 
Only  dte  third  body  efficiency  for  Ar  is  precisely  known.  As  seen  on  the  reacUon  pathway  diagram  in 
Figure  8,  lire  rates  of  reactions  for  H2O  aird  NjO  as  collision  parmers  are  most  significam,  aiKl 
NnO  tiavc  llK  largest  effect  on  ilte  NjO  oonocniratiocr  as  indreatsd  in  the  ,N,0  scstsiiivity  plot  for  tits 
N2O  M  reaction  tor  different  coUiders.  shown  in  Figure  1 1.  The  rcactzuit  as  a  collider  does 
appear  im[)ortaiii  possibly  due  to  the  low-pressure  c«>aitions  whiclr  auisitcc  the  di {fusion  of  light  species, 
thus  reducing  the  concentration  of  Kj  trear  the  bumcr  surface.  However,  at  higher  {Ncmires  diffusion  is 
not  so  large  an  effect  and  as  a  collider  may  becoore  more  imponani.  Also,  tire  efficiency  for  Hj  as 
a  colltder  is  assumed  to  be  small,  as  is  frequently  observed  for  many  reactions.  This  assumption 
cocitnbotcs  to  the  prediction  of  low  seitsitivity  for  the  reaaion. 


4.3J  Conv'ersion  of  NjO  to  NO.  Formation  of  a  few  pcrccni  of  NO  has  an  traporuttt  inhibiting 
effect  OR  the  overall  rare  of  ownbustion  because  less  heat  is  released  w=hcn  it  is  formed  compared  to  the 
more  typical  products.  Nj  and  KjO.  Presented  in  Table  2  is  a  w^mporison  of  the  measured  butrxsd  gas 
composition  for  the  present  flame  at  25  mm  above  the  burner  surface  with  that  calculated  with  tfjc  flame 
code.  Also  pre-sented  in  3'able  2  arc  Ue  species  concentratiotts  obtairred  from  a  NASA-Lewis  chemical 


NORMALIZED  SENSITIVITIES 


Table  2.  Experimental  and  Calculated  Major  Species  Concentrations  in  the  Post  Rame  Region®: 
Units  are  in  mol/cm^  (x  10“®) 


Species 

Experiment 

Premix  Model 

Nssa-Lewis 

..... ,  - - 

(Equilibrium) 

. . . . . . 

li. 

0.45 

0.546 

0,154 

NjO 

0.00 

0.0  J  6 

0.000 

N2 

5.4 

5.20 

5.53 

HjO 

4.4 

4.81 

5.35 

NO 

0.45 

0.450 

o.on 

Oj 

0.032 

0.038 

0.060 

*  Height  tbove  btancf  i5  ntjn  *n<l  tcsapartCtfe  -2,000  K, 


eq>iilibriuai  code  (Svdtla  and  McBride  1973)  calculation  using  a  temperature  of  2,0(X)  K.  Ttiis 
temperature  corresponds  to  a  height  of  about  25  ram  in  the  fkme.  Noiitx  that  the  measured  concentrations 
and  those  calculated  wiiii  the  tlamc  code  agree  ^masonably  well.  However,  both  tire  measured  and 
predicted  concemratiorts  of  aixl  NO  arc  considerably  larger  ihatt  tlic  equilibrium  concenunaiic'ns. 
Careful  inspection  of  the  tabic  reveals  that  H^O  and  Oj  arc  predicted  by  the  flame  model  to  be  slightly 
lower  than  equilibrium.  TlUs  result  is  best  seen  by  comparing  tire  calculated  concentrations  since  tlw; 
diffcrcTKie  is  much  smaller  than  the  precision  of  die  measytcmenis.  Using  iIk  predicfal  cottcenuatiom 
from  the  flame  code  and  the  measured  temperature  of  7.,0iM)  K  input  parameters  to  the  equilibrium  code, 
the  calculated  adiabatic  flame  temperature  of  this  mixture  is  found  to  bs  ap«.>roxin)stely  2,200  K,  21X)  K 
higlicr.  As  further  evidence  of  the  irapottaitcc  of  Ute  NO  formation,  flame  speeds  were  calculated  for  fljc 
free  flame  corresporrding  to  our  conditions.  (Note,  however,  cherc  arc  no  heat  losses  in  the  energy 
equation.)  For  the  mechanism  givat  in  TaWc  1.  dte  calculated  Rame  speed  is  179.0  cm/S.  When  all  of 
the  NO  tcaaions  are  removed,  trie  calculated  flame  speed  is  22S.2  ctn/s.  Thus.  c««  socs  that  the  fomtaiion 
''  *■  NO  is  not  only  of  signiflcance  irr  poDutani  formatjcn.  but  can  also  Isc  very  iraportsnt  in  deteiriuning 
the  gross  behavior  of  NRQ  oxidired  combustion. 

A  sensitivity  plot  for  [10  obtained  using  the  ntc'dranism  in  Tabic  1  is  sisown  in  Hgure  12.  The 
reaction  NjO  K  =  NH  +  NO  <-R20)  is  highly  .sensitive  it>  Ute  burned  gas  regioti  and  is  the  major  sourec 
of  NO  ptoduaion  at  the  presem  cursditions.  Until  recetidy,  the  rate  expression  was  tsot  well  csiablisJicd 
frem  the  few  high-tempcreiuie  roeasenrments  available  since  all  of  fliesc  involved  comr>lex  mechanism 
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fits.  Consequently,  a  wide  range  of  rate  expressions  has  been  used  to  model  N2O  combustion  (Cattolica, 
Smooke,  and  ESean  1982;  Coffee  1986;  Miller  and  Bowman  1989;  Miller  ei  al.  1983).  Two  of  the  rate 
expressions  used  by  Miller  and  coworkers  (Miller  and  Bowman  1989;  Miller  et  al.  19S3)  are  considerably 
larger  than  the  others  in  ute  temperature  range  of  interest.  Initial  attempts  to  model  tlie  present  experiment 
using  the  expression  of  Miller  and  Bowman  (1989)  resulted  in  an  underprediction  of  the  NO  mole  fraction 
and  an  overprediction  of  O2  in  the  burned  gases  by  approximately  25%  each.  The  computed  O2 
concentration  in  die  burned  gas  region  has  a  strong  negafive  sensitivity  for  R20  (see  Figure  10)  because 
of  NH  production  followed  by  subsequent  reactions  R22  and  R23  and  tte  sequence  -R27  followed  by 
“R37.  Using  the  earlier  expression  of  Miller  et  al.  (1983)  for  R20,  which  is  slightly  larger  at  high 
temperrtures,  brought  the  production  of  NO  and  Oj  to  within  5%  of  the  experimental  values.  A  recent 
transition  state  theory  calculation  by  Mider  and  Melius  (1992)  has  led  to  a  theoretical  expression  which 
is  similar  to  that  of  Miller  et  al.  (1983)  and  is  our  recommendation.  Also,  a  recent  experimentai  and 
modeling  study  oy  Martin  and  Brown  (1990a,  1990b)  on  CHyair  and  H2/02/Ar  flames  doped  witii  NjO, 
NO,  or  NH3  strongly  supports  our  recommend^oa  In  their  work,  detailed  rate  ^  sensitivity  analyses 
were  performed  for  the  N^O  in  the  two  flames  containing  NH3.  It  was  shown  that  the  NjO  in  these 
flames  is  formed  almost  exclusively  by  R20  and  destroyed  primarily  by  R15  (P19  was  also  considered). 
By  examining  the  sensitivity  coefficients  presented  in  that  work  for  N2O,  one  can  show  that  if  a  much 
smaller  rate  coefficient,  say  a  factor  of  5-10,  were  to  be  used  for  R20,  the  N2O  concentraticsis  would  be 
severely  undeipredicted.  Thns,  it  appears  the  rate  of  -R20  is  much  laiger  than  that  used  by  Cattolica, 
Smooke,  and  Dean  (1982)  and  Coffee  (1986). 

The  reaction  of  NH  with  NO  has  been  the  subject  of  a  number  of  theoretical  and  experimental 
investigations  (Gonden,  Mulac,  and  Nangia  1971;  Hansen  et  ai.  1976;  Harrison,  Whyte,  and  Phillips  1986; 
Yamasaki  et  al.  1991;  Vandooren  et  al.  1991;  Mertens  et  al.  1991;  Yokoyama,  Sakane,  and  Fueno  1991; 
DuraiU  and  Rohlfing  1993;  Michaud,  Westmoreland,  and  Feitelburg,  in  press;  Bozzclli  and  Dean  1993; 
Fueno,  Fukuda,  and  Yokoyama  1988;  Harrison  and  Maclagan  1990),  most  of  these  very  recent,  which 
center  on  the  rate  of  reaction  at  various  temperatures  aid  the  channels  for  product  formation.  There  is 
remarkably  good  agreement  among  the  low-tempemture  determinations  (Gordon,  Mulac,  and  Nangia  1971; 
Hansen  et  ai.  1976;  Harrison,  Whyte,  and  Riillips  1986;  Yamasaki  et  al.  1991),  so  that  the  overall  rate 
coefficient  is  near  collisional  with  a  value  of  approximately  3  x  10^^  cm^/mol-s.  Tliis  result  is  invariant 
at  a  pressure  range  from  at  least  300  ton  (Gordon,  Mulac,  and  Nangia  1971;  Hansen  et  al.  1976;  Harrison, 
Whyte,  and  .Phillips  1986;  Yamasaki  et  ai.  1991),  and  using  a  wide  array  of  carrier  gases  (cf.  Harrison, 
Whyte,  and  Pt/illips  1986).  At  higher  temperatures  the  agreement  between  the  experiments  is  fair,  but  not 
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as  good  (Vandoorcn  ct  al.  1991;  Mcricns  ct  ai.  1991;  Yokoyama,  Sakane,  and  Fueno  1991).  Of 
experiments  wherein  the  branching  ratios  were  measured,  three  favor  ihe  NjO  +  H  channel  (Mettens  et  al. 
1991;  Yokoyama,  Sakane,  and  Fueno  1991;  Durant  and  Rciilfing  1993)  in  agreement  with  our  work,  and 
two  favor  the  N.2  +  OH  chanr  d  (Yamasaki  et  al.  1991;  Vandooren  et  al.  2991).  The  work  of  Vandoren 
et  al.  (1991)  was  an  interpretation  of  an  NH3  flame  experimmit  wherein  the  results  are  highly  dej^ndent 
upon  absolute  concentrations  of  several  of  the  radical  species  which  could  not  be  measured  better  than  to 
within  a  factor  of  two.  The  resulting  rate  coefficients  must  have  very  large  error  limits,  quite  possibly 
large  enough  to  be  compatible  with  N2O  actually  being  ifse  preferred  product  'Hre  total  of  the  two 
channels  is,  however,  in  fair  agreement  with  ofher  high-temperature  measukc.ments.  Y amssaki  et  al.  (1991) 
concluded  that  only  the  N2  +  OH  products  are  fonned.  It  is  not  clear  why  this  disagrees  so  strongly  with 
the  other  results.  In  addition  to  these  experimental  results,  there  have  been  several  theoretical  calculations 
(Miller  and  Melius  1992;  Durant  and  Rohlfing  1993;  Michaud,  Westmoreland,  and  Feitelburg,  in  press; 
Bozzelli  1993),  all  of  which  conclude  tha  die  N^O  channel  is  stiujgly  preferred.  The  main  reason  for 
this  involves  the  energetics  of  the  reaction’s  potential  surfaces.  All  of  the  surfac'  studies  on  the  reaction 
(Miller  and  Melius  1992;  Durant  and  Rohlfing  1993;  Fueno,  Fukuda,  and  Yokoyama  198)i;  Harrison  and 
Maclagan  1990;  Walch  1993)  agree,  qualitatively,  diat  the  reaction  starting  from  NH  +  NO  proceeds 
through  a  cis  HNNO  intermediate  in  a  state  which  connects  adiabadcaily  with  the  two  product 
channels.  Tlie  barrier  to  H  atom  transfer,  leading  to  NNOH  which  rapidly  dissociates  to  N2  +  OH,  is 
approximately  6  kcai/wol  higher  than  that  for  simple  H  atom  elimination.  The  reaction  therefore  strongly 
favors  the  H  +  N2O  channel.  Walch  (1993)  has  carefully  examined  the  NH  +  NO  entrance  and  found  that 
there  are  barriers  of  a  few  kcal/mol  in  tne  ^A'  surface,  but  none  in  the  other  symmetry  allowed  surface, 
^A".  The  surface  leads  to  an  excited  state  of  HIMNO  which  could  stabilize,  but  which  cannot  further 
react  to  products.  In  contrast,  the  low-temperature  study  of  Harrison,  Whyte,  and  Phillips  (1986)  indicates 
that  there  is  no  appreciable  barrier  to  the  overall  reaction.  The  observed  near-collisional  rate  at  low 
temperature  also  strongly  supports  this  conclusion.  However,  the  fact  that  the  overall  rate  exldbits  no 
pressure  dependence  shows  that  .stabilized  liNNO  cannot  be  the  major  product.  If  Walch  (1993)  is  correct 
concerning  (he  entrance  channel  barkers,  these  results  suggest  that  the  favored  reaction  patli  involves 
approach  of  NH  and  NO  on  one  of  the  ^A"  surfaces  followed  by  a  crossing  to  the  cis  ^A'  surface  on 
which  the  final  products  can  be  reached.  Corroborative  studies  of  these  entrance  channel  characterisdes 
and  die  crossing  probability  would  be  nost  interesting. 


The  NO  sensitivity  plot  in  Figure  12  indicates  that  several  reactions  involving  the  NH  and  HNO 
radicals  and  N  atoms  have  a  more  modest  effect  Oii  NO  formation.  This  result  may  be  understood  by 
considering  the  reaction  pathway  diagram  in  Figure  8.  Some  of  the  NH  radicals  formed  by  -R20  are 
converted  to  N  atoms  by  mactions  with  H  or  OH.  The  N  atoms  primarily  react  vrith  NO  to  form  N2, 
although  some  form  additional  NO  by  a  reaction  with  OH  (~R29).  The  mechanism  has  a  small  portion 
of  the  NO  fonrned  via  the  padiway  NH  ^  HNO  ^  NO  involving  two  OH  reactions.  This  pathway  for  NO 
formation  is  probably  the  most  poorly  established  in  ilie  mechanism  tecause  the  rate  constant  for 
NH  +  OH  =  HNO  4-  H  (R24)  may  be  traced  to  an  estimate  (Miller  et  al.  1983). 

llie  reaction  N2O  +  NH  =  HNO  +  N2  (R41)  was  also  considered  in  the  mechanism.  Its  rate  constant 
is  not  well  established  and  the  expression  used,  k  =  2.0  x  I0^^e~^®’®®^^cm^/mol-s,  was  taken  from  an 
estimate  tabulated  in  Hanson  and  Salimian  (1985).  After  consideration,  it  was  excluded  from  the  leactitm 
set  since  no  changes  in  the  major  species  profiles  were  observed  with  its  inclusion.  There  is,  however, 
a  slight  increase  in  the  NO  concentration  because  of  the  conversion  of  HNO  to  NO  '"nd  a  subsequent 
decrease  in  N  atoms  due  to  tiseir  reaction  with  NO.  Reaction  R41  has  the  largest  effect  on  the  HNO 
concentration  which  doubles  when  the  reaction  is  included  in  the  mechanism.  The  effects  of  R41  on  the 
HNO  and  N  profiles  arc  shown  in  Figure  13.  The  importance  of  R41  in  the  mechanism  is  difficult  to 
assess  since  the  HNO  and  N  radicals  could  not  be  measured  in  the  present  sti  dy  and  none  of  the  measured 
.species  has  much  sensitivity  to  the  teaction.  Therefore,  the  relative  importance  of  R41  deserves  future 
consideration. 

5.  CONCLUSION 

A  combined  expen  rental  and  detailed  chemical  modeling  study  of  a  stoichiometric  H2/N20/Ar  flame 
has  been  presented.  TIvj  experiments  were  performed  at  low  pressure  with  MB/MS  ai>d  LIF  diagnostic 
techniquei,.  The  low-pressure  environment  enables  the  profiles  to  be  examined  at  higher  spatial  resolution 
than  has  appeared  in  most  of  the  prior  flame  studies  on  this  chemical  system.  In  addition,  a  larger  number 
of  species  has  been  measured  than  in  previous  studies. 

A  cnemical  mechanism  for  the  system  was  developed  by  compiling  a  comprehensive  mechanism  and 
thi  n  remc  ing  unimportant  species  and  most  of  the  unimportant  reactions  tiirough  sensitivity  aruilysis. 
Festinent  ll^ratuie  was  reviewed  in  order  to  select  the  best  rate  coefficients  tor  the  important  reactions. 
In  general,  tlte  agieetnent  between  the  model  and  experiment  is  acceptable.  The  chemistry  governing  the 
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present  flame  conditions  has  been  discussed  in  detail.  The  main  qualitative  features  of  the  mechanism 
have  not  changed  a  great  deal  from  previous  recommendations.  However,  several  important  new  findings 
liave  been  made.  First,  the  possible  role  of  the  reaction  N2O  +  OH  =  HO2  +  N2  has  been  investigated. 

present  results  arc  best  modeled  if  the  rate  coefficient  used  is  very  small;  however,  a  modest  rate 
coefficient  may  be  acceptable.  If  it  is  fast  enough,  this  reaction  could  have  important  effects  on  die 
overall  combustion  rate  under  pertinent  conditions.  Second,  the  importance  of  various  colliders  in  the 
reaction  NjO  +  M=N2  +  0  +  M  was  also  considered.  Many  previous  studies  have  established  that  this 
reaction  is  a  major  initiation  step  in  the  system.  Under  the  present  conditions,  the  most  important  colliders 
appear  to  be  N2O  and  H2O.  However,  this  result  depends  in  part  on  the  supposititm  that  the  estimated 
efficiency  factor  for  H2  is  considerably  smaller  than  that  of  N2O  and  H2O.  Third,  the  present  work 
corroborates  the  results  of  other  recent  studies  which  indicate  that  the  rate  coefficient  of  the  reaction 
NH  +  NO  =  N2O  +  H  should  be  larger  at  high  temperature  than  indicated  in  previous  studies. 


lNTENn(»iALLY  LEFT  BLANK. 
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